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Abstract 



We investigate the polarization fluctuations caused by quantum noise in quan- 
tum well vertical cavity surface emitting lasers (VCSELs). Langevin equations 
are derived on the basis of a generalized rate equation model in which the in- 
fluence of competing gain-loss and frequency anisotropics is included. This 
reveals how the anisotropics and the quantum well confinement effects shape 
the correlations and the magnitude of fluctuations in ellipticity and in po- 
larization direction. According to our results all parameters used in the rate 
equations may be obtained experimentally from precise time resolved mea- 
surements of the intensity and polarization fluctuations in the emitted laser 
light. To clarify the effects of anisotropics and of quantum well confinement 
on the laser process in VCSELs we therefore propose time resolved measure- 
ments of the polarization fluctuations in the laser light. In particular, such 
measurements allow to distinguish the effects of frequency anisotropy and of 
gain-loss anisotropy and would provide data on the spin relaxation rate in 
the quantum well structure during cw operation as well as representing a new 
way of experimentally determinig the linewidth enhancement factor a. 

PACS numbers: 42.55.-f, 
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I. INTRODUCTION 



Vertical Cavity Surface Emitting Lasers (VCSEL) have recently attracted intensive ex- 
perimental and theoretical effort. One of the major advantages VCSELs have over conven- 
tional semiconductor lasers is the highly symmetric geometry around the axis of laser light 
emission. The polarization of light emitted from VCSELs is therefore not determined by the 
massive anisotropy of the device architecture, as is the case in edge emitting lasers. Instead, 
due to the transverse symmetry of the cavity and the active region the polarization is highly 
sensitive to more subtle effects, such as small anisotropies in the crystal structure, strain or 
optical anisotropies in the mirrors [Jlj-Q • 

Experiments have already provided a large range of results on polarization as a function 
of pumping current @J§]. The variety of results indicates the sensitive dependence of po- 
larization on very small effects. However, the connection between the observed stability or 
bistability of linear polarization and the theoretical explanations offered are still no more 
than tentative. 

In 1995, San Miguel, Feng and Moloney introduced a rate equation model for quantum 
well VCSELs, based on the observation, that the electron-hole pairs in the bands closest to 
the band gap can be separated into electron-hole pairs emitting only right circular polarized 
light and electron-hole pairs emitting only left circular polarized light J7[. Effectively, this 
corresponds to two independent reservoirs of electron hole pairs, each characterized by its 
own electron-hole pair density, coupled only by spin relaxation processes. In the following, 
this assumption will be refered to as the split density model. 

The first major result derived from this model has been the proposal that polarization 
stability may be a consequence of a frequency anisotropy due to birefringence in the optical 
cavity ||. As opposed to the more straightforward assumption of a gain- loss anisotropy, 
the frequency anisotropy may give rise to polarization switching as the injection current is 
increased. However, since polarization switching may also be induced by temperature related 
effects 0, the fact that polarization switching is experimentally observed is not sufficient to 
verify the assumption of dominating frequency anisotropies. Although further work based 
on the split density model of San Miguel and coworkers has been recently forthcoming , 
the predictions made by the model have so far been insufficient to truly rule out the presence 
of alternative mechanisms. In particular, no attempts have yet been made to determine the 
anisotropies and the time scales involved in the laser process directly from the laser emission 
during cw operation. 

In order to demonstrate that such a direct test of the model is indeed possible, we inves- 
tigate the polarization fluctuations predicted by the Langevin equations derived from the 
split density model. One of the advantages of this approach is that polarization fluctuations 
can be measured at a fixed injection current during cw operation. Deriving the anisotropies 
from fluctuation measurements provides not only a tool for testing the theory, but can also 
provide information about temperature dependent changes in the anisotropies. In this way, 
our theory may help to resolve questions such as those raised in and |J. 

In section [TI]of this paper, we formulate the full rate equations of the split density model, 
using the normalized Stokes parameters to describe light field polarization. In section JT|, we 
examine the most likely case of a frequency anisotropy and a gain-loss anisotropy along the 
same or two orthogonal crystal axes. The linearized Langevin equations for the stationary 
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point are derived. In section [Ty|, the solution of the Langevin equation is given, using 
reasonable approximations with regard to the time scales involved. In section 0, the results 
are discussed and possible experiments are proposed. Finally, in section [VT|, conclusions and 
an outlook are presented. 



II. FORMULATION OF THE RATE EQUATIONS FOR THE SPLIT DENSITY 

MODEL 

A. The split density model 

If only the lowest lying conduction bands and the highest lying valence bands contribute 
to the laser process in a quantum well structure, the preservation of angular momentum 
around the axis perpendicular to the quantum well and the fact that photons which are 
emitted along this axis have angular momenta of ±1 (corresponding to either left or right 
circular polarization) limits the number of possible emission processes to two separate tran- 
sitions. Conduction band electrons with a spin of +1/2 around the axis perpendicular to 
the well can only recombine with heavy holes of -3/2 angular momentum around this axis, 
emitting a photon with an angular momentum of -1. (Note, that a hole with an angular 
momentum of -3/2 corresponds to an empty electronic level with an angular momentum 
of +3/2.) Correspondingly, spin -1/2 electrons only recombine with holes of +3/2 angular 
momentum, emitting photons with an angular momentum of +1. 

This means that effectively two completely distinct pools of electron-hole pairs exist. In 
one pool each electron-hole pair has an angular momentum of -1 and interacts only with 
right circular polarized light. The electron-hole pairs in the other pool, having an angular 
momentum of +1 per pair, interact only with left circular polarized light. Rate equations 
based on this assumption have originally been formulated by San Miguel and coworkers in 
0. In the following we adopt the notion of the split density model and considerably extend 
the analysis towards a generalized representation able to accomodate in the model arbitrary 
types of anisotropies. 



B. Parameters 

As in the paper by San Miguel 0, we use the variables D for the total electron-hole 
density and d for the difference between the densities of electron- hole pairs with +1 and 
with -1 angular momentum. We will use a somewhat different approach to describe the light 
field in the cavity, however, to emphasize the difference between light field intensity and light 
field polarization. The total number of photons in the cavity (regardless of polarization) is 
denoted by n. The polarization is then described by a three dimensional vector of unit length, 
P, which defines the point on the Poincare sphere corresponding to the present polarization. 
The components of P are the normalized Stokes parameters. If E± are the complex electric 
field amplitudes of right and left circular light, the normalized Stokes parameters are given 
by 

P ~ g±g= + E - E+ (la) 
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The spatial directions of the vector P correspond to pairs of orthogonal polarizations. The 
component of P along any given direction is equal to the intensity difference between the 
two polarizations relative to the total intensity. P3 is the intensity difference between left 
and right circular polarization, P\ is the intensity difference between x and y polarization 
and P2 is the intensity difference between light polarized along the (1,1) direction and light 
polarized along the (1,-1) direction. In the plane of linear polarization, a rotation of P by 
180 degrees therefore corresponds to a 90 degree change in the direction of polarization. 
While this may not seem a very intuitive picture at first, one of the advantages of describing 
the light field in terms of P is that it directly refers to light field intensities as they are 
observed in experiment. 



C. Timescales and anisotropics 

The laser process in the split density model for quantum well VCSELs is characterized 
by four timescales and three anisotropics. From slowest to fastest, the timescales are defined 
by 

1. The rate of spontaneous emission into the laser mode 2w, usually around 10 6 s _1 to 
10 7 s _1 , although this can be very much a function of cavity design [10,11]. 



The rate of spontaneous carrier decay 7, usually around 10 9 s _1 to 10 10 s _1 , which also 
depends strongly on the cavity. In fact, the reduction of this decay into non-laser 
modes is at the heart of a proposal to achieve thresholdless lasing by optimizing the 
optical cavity |TT . 



The rate of spin relaxation 7^. This rate is an unknown quantity, especially in the 
high carrier density regime crucial to the laser structures under investigation. In this 
model, it appears only as the sum with the spontaneous carrier decay which we define 
as T = 7 S + 7. Experiments and theory indicate, however, that the order of magnitude 



will probably be similar to that of 7 [12|- [13]. 



4. The rate of photon emission from the cavity 2k. 

The factor of 2 in the rate of spontaneous emission is a result of using photon number and 
intensity variables instead of the field amplitudes to describe the light in the cavity. Also, 
note that both w and k are averages over anisotropic properties of the VCSEL. 

The three possible anisotropies in gain, loss and frequency are characterized by both 
their magnitude and by their geometrical orientation. Similar to the definition of the Stokes 
vector P, it is possible to represent the anisotropies by vectors. The direction of this vector 
corresponds to the direction of P for which the extremal gain/loss/frequency values are 
obtained. The length corresponds to the difference between the extremal values. The 
anisotropy vectors are defined as 
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1. gain anisotropy g, such that the rate of spontaneous emission is given by 2w(l + Pg). 

2. loss anisotropy 1, such that the rate of photon emission from the cavity is given by 
2k(1 + P1). 

3. frequency anisotropy ft, such that the length of Q is equal to 1/2 the frequency differ- 
ence between the modes of orthogonal polarization. 



We can now formulate the rate equations for any arbitrary set of anisotropies using the 
parameters defined above. 



e 3 indicates the unit vector in the direction of the 3rd component of the stokes vector. 

Dq is the injection current in units of 7 and a is the linewidth enhancement factor, which 
describes a shift in frequency due to the electron-hole density in the quantum well. In the 
presence of two separate reservoirs of electron- hole pairs (split density model), where the 
carriers from one reservoir interact only with one type of circular polarization, the frequency 
shift causes an effective birefringence ||. 

To understand the terms in the equation, it is important to keep in mind that each time 
n or w appear in the equations, they are modified by the anisotropy factors of (1 + PI) and 
(1 + Pg), respectively. In the model this is the only influence of the anisotropies on the 
dynamics of D and n if d — 0. If d ^ 0, the last terms in (2a) and (2b) increase or decrease 
the rate of stimulated emission, depending on whether the light field in the cavity interacts 
more strongly or more weakly with the available carriers. The last term in (2c) describes 
the hole burning effect of P3 7^ 0. 

The dynamics of the Stokes parameters P is directly determined by the anisotropies. 
The frequency anisotropy induces a rotation around the axis defined by O. The gain-loss 
anisotropy draws P towards one of the poles along the axis defined by g and 1, at a rate 
proportional to the component of g or 1 orthogonal to P. This is illustrated in Fig. |l|. 

The effect of the split density model on the polarization dynamics can be understood in 
terms of gain and frequency anisotropies introduced by the density difference d: 



D. Rate equations 




(2d) 
(2e) 



(2b) 



(2a) 



(2c) 



d 

Seff = g + -^e 3 
n eff = ft + w(l + Pg)ade 3 . 



(3) 
(4) 
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In passing we note that another consequence of the split density model is the fact that the 
factor of 2 in the rate of spontaneous emission is absent. This is due to the two parallel laser 
processes having an induced emission rate of 2w(D/2)(n/2) for d = and P 3 = 0. The sum 
of the two emission rates is therefore wDn, as given above. 

III. LANGEVIN EQUATIONS FOR CW OPERATION AT A STABLE LINEAR 

POLARIZATION 

A. Anisotropies along the [110] and [110] crystal axes 

Although a large number of nonlinear effects and dynamical properties can be described 
by chosing special combinations of anisotropies || , we will concentrate on the simple case 
of cw operation at a stable linear polarization. For many VCSELs, this seems to be the 
natural state of affairs, even when no artificial anisotropies were created during the growth 
of the device. It has been found that a large number of VCSELs emit light polarized along 
the [110] or [110] crystal axes - @ . This must obviously be the result of unintentionally 
introduced anisotropies. Possible reasons for this are the slight tilt of the growth axis often 
used in MOCVD [|] or the tendency of strain induced changes in the optical properties 
to produce anisotropies along these axes || . The latter effect can be visualized quite 
nicely. If one looks at a [001] surface of a semiconductor lattice, the projections of the bonds 
appear along the [110] and the [110] directions. Light polarized along one of these directions 
will mainly interact with the electrons in the bonds along this direction. Consequently, 
the largest part of the optical anisotropy induced by stress is caused by the difference in 
compression between the [110] and the [110] bond directions. 

Given this preference for two orthogonal linear polarizations observed in many VCSELs, 
we choose g = ge 1: 1 = /ex and O = Qe-i. The rate equations are now 

r (D-D )-w(l + P ig )dnP 3 (5a) 
(1 + P x l)n + w(l + Pig)dnP 3 (5b) 
d-w(l + P x g)DnP 3 (5c) 

2k(1 + Pil)l){Pl - 1) - w(l + P 1 g)d(P 3 P 1 - aP 2 ) (5d) 
2«(1 + Pil)l)PiP 2 - w(l + P x g)d{P 3 P 2 + aPi) - fiP 3 (5e) 
2«(1 + Pil)l)PiP 3 - w(l + Pig)d{Pl - 1) + QP 2 . (5f) 

B. Linearization around the stationary solution 

The stationary solution of these equations is given by Pi = 1, P 2 = P 3 — d — and 
D = 2k(1 + l)/w(l + g). The photon number in the cavity n s is a linear function of D , 
specifically 



-D = -w(l + P 1 g)Dn-~ l 
dt 

d 

—n = w(l + Pig)Dn — 2k 

(JjL 

-f-d = —w(l + P\g)dn — V, 
dt 

j t P l = -(w(l+P 1 g)Dg- 
j t P 2 = -( w (l + P ig )Dg- 
±P 3 = -(w(l+P ig )Dg- 
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Close to this stationary point the laser relaxation dynamics can be linearized. Deviations 
from the stationary point are described by five coupled dynamical variables. In the case un- 
der consideration the dynamics can be seperated into two mutually decoupled sub-systems. 
One sub-system describes the dynamics of the total photon number fluctuations 5n = n — n s 
coupled to the total density fluctuations 5D = D — 2k(1 + 1) / w(l + g) . The other sub-system 
describes the coupling of the density difference d to the polarization parameters P<i and P3. 

To obtain the linearized Langevin equation, we further add the noise term f(t) to the 
dynamics. This time dependent five dimensional vector incorporates all external noise and 
signals influencing the laser. In the case under consideration, this will be the vacuum 
fluctuations of the light field entering the cavity. A discussion of the statistical properties 
of this noise term will be given in the next subsection. 

At this point, it is convenient to introduce five new parameters for a more compact 
formulation: 

1. The injection current in units of the threshold current is defined as x, which means 
that the stationary photon number n s is replaced by 

x = ^±^n s + 1. (7) 

7 

2. The gain-loss anisotropy is combined into a single variable p equal to the difference of 
the two, scaled with the ratio of cavity loss rate and spontaneous carrier decay. This 
dimensionless quantity should be of the order of unity to be effective, corresponding 
to a relative anisotropy of about 0.1 percent for typical timescales: 

25(1 + (8) 

7 

3. The frequency anisotropy is scaled in terms of the spontaneous carrier decay rate 
divided by a, as it only stabilizes the polarization in conjunction with a [|J: 

9 = ^. (9) 

7 

4. The ratio of spin relaxation rate and spontaneous carrier decay rate is written as r, 
such that 

r , , 

r = 1. (10) 

7 

5. The relaxation oscillation frequency v is defined as a function of the injection current 
x: 



V = ^2/c(l + lh(x - 1) (11) 
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Using these parameters, the Langevin equation for the quantum well VCSEL is 
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or at a rate of v 



(12) 



it is realistic to assume that 7 <C k, the spontaneous carrier decay rate 7 will usually 
be significantly smaller than the relaxation oscillation frequency v = ^2k{1 + Vyy{x — 1). 
Note that this condition breaks down very close to threshold (x = 1). However, as close 
to threshold the relative polarization fluctuations approach infinity there is no polarization 
stability for x very close to one. 



C. Quantum noise 

Although the formulation given above can be applied also to problems of externally 
injected fields and similar linear response problems, we will now define the noise term f(t) 
as the electro-magnetic field noise entering the cavity from the vacuum. 

The definition of a noise term for the light field in a cavity is a standard procedure 
in quantum optics. In the case of two modes of orthogonal polarization, the extension is 
straightforward. The two time correlation functions which are non-zero are 



(fsn(t)fsn(t + r)) =4k(1 + ZH<S(t) 
4«(l + Z)7(a;-i; 



w(l + g) 



S(t) (13a) 



7 (x-l) V) V ' 

(f P3 (t)fp 3 (t + r)> = _ S(t). (13c) 

To understand the derivation of these terms, one may either think in terms of a classical 
field, considering the vacuum modes and the dipole densities associated with the carrier 
density pools as field modes with random quantum fluctuations acting as external forces on 
the field in the cavity, or one may apply particle picture reasoning: since only whole photons 
are emitted into and out of the cavity, there is a stochastic process involved which gives rise 
to shot noise. It is one of the fascinating properties of quantum mechanics, that the same 
noise terms result from two seemingly different pictures. 

Another noise source for the laser process is the noise in the carrier injection. However, 
this noise is much weaker than the light field noise if •yx k, which is a natural assumption 
for real devices. The physical reason for the relative smallness of shot noise from injection is 
that, if the photon emission rate from the cavity is much faster than the total carrier decay 
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rate, then the number of carriers present in the active region is far greater than the number 
of photons in the cavity. Consequently, the relative statistical fluctuations are much smaller 
in the carrier subsystem than in the light field. 

Note that this approximation is only valid for short timescales. On timescales much 
longer than the I/7, energy conservation requires the fluctuations described here to cancel. 
This means that the low frequency noise of the power spectrum is very sensitive to the weak 
fluctuations in the carrier injection [I5|. 



IV. SOLUTION OF THE LANGEVIN EQUATION 



A. The Greensfunction solution near the stationary point 



In accordance with the fluctuation-dissipation theorem, the noise in the light emitted 
from the laser is approximately given by a linear response to the quantum noise entering 
the laser cavity. The fluctuations can therefore be calculated from the linear response of the 
VCSEL. The five dimensional Greensfunction can be obtained by determining the Eigen- 
values and Eigenvectors of the non-symmetric five by five matrix describing the linearized 
dynamics. Since the Eigenvectors are non-orthogonal, both left and right Eigenvectors need 
to be determined. 

The problem separates into one two dimensional problem and one three dimensional 
problem. Therefore, an exact analytical solution is possible. However, to understand the 
physical significance and to single out the experimentally relevant case, it is useful to apply 
the realistic assumption that, since 7a; <C k, 7 <C u, except for x very close to 1. Also, 
the anisotropies p and 9 should not be much greater than 1. Using these assumptions, the 
Eigenvalues Aj and Eigenvectors a, and bj are 



A1/2 = ~\ix ± 



w 



ai/2 



1 

V2 



7(*-l) 



10 



h/2 



1 

V2 



7(3-1) 
1 









(14) 



A 3 = -j{p + 6) 



a 3 = ( 1 -a ) b 3 






1 



(15) 



A4/5 = -f 7(2 + r + p — 9)±iv 
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^4/5 = ^(0 T ^ 1 



>4/5 
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±i 



w(l+g) 

a 



(16) 



Using this set of vectors, any external perturbation of field variables or carrier densities 
can be decomposed into the Eigenvectors with the corresponding exponentially decaying 
and oscillating Greensfunction. The time integral over the perturbation then gives the 
linear response of the laser dynamics. This procedure can now be applied to quantum noise. 



B. Application of the Greensfunction to quantum noise 

The left Eigenvectors aj are used to decompose the noise into contributions associated 
with the corresponding Eigenvalues. If the noise terms are represented in the form of a five 
by five diffusion matrix N [JTJJ], the decomposition can be accomplished by calculating the 
matrix elements iVjj = ajNaj. The fluctuation matrix F(r) can then be expressed as a sum 
over the dyadic products of the right Eigenvectors b«. 

FW = E^n;^®bj (i7) 

As the Eigenvectors are not orthogonal, there are contributions for i ^ j with complex 
— \i — A*. Since we apply the approximation that 7 z/, however, the absolute value 
of these terms is much smaller than those with i = j and we can obtain the light field 
fluctuations by summing over the five terms with i = j only. This effectively means that the 
Eigenvectors fluctuate independently. The fluctuation matrix is then given by a sum over 
projection operators which must be symmetric. 



C. Light field fluctuations 

The light field fluctuations of polarized light have three degrees of freedom, given by the 
total intensity n, the ellipticity P3 and the direction of linear polarization P2. Including 
possible correlations, these fluctuations are described by a symmetric three by three matrix. 
However, since the dynamics of the total total intensity is decoupled from the polarization 
dynamics, there is no correlation between the total intensity and the polarization. Therefore, 
four fluctuation terms are sufficient to completely describe the fluctuations in the laser light 
during cw operation: 

(5n(t)5n(t + r)> = ^ + - ~ V ^cosfrr) (18) 
v w v " w (1 +g) x v ' v ' 

(8n(t)8n(t + r)) A /2 

- — — — ^ — = — -e 1 ■ ' cos(ut) (19) 

n 2 s x(x-l) K J K ' 
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(P 3 (t)P 2 (t + t)> = a(P 3 (t)P 3 (t + t)) 

a ^ -7(s+r+p-e)r/2 



(x — l)(x + r + p 



e-^ x+r+p -°> T ' 2 cos{vT) (21) 



{PMP * + T)> = (i - l)(x + r + p- e f " fl ' + ' + ^"' Tl2ms (' /T ) 
The factor A is a measure of the overall magnitude of noise. 

a = Mi±jMi±g) (23) 

7 2 

Fig. shows the noise distribution at r = for a realistic choice of parameters. Fig. |3| shows 
the polarization fluctuations as a function of r. 

There is a correlation between the two fluctuations, mediated by a, which is a typical 
feature of the two density model. The a factor converts the density difference fluctuations 
d into frequency difference fluctuations between left and right circular polarization. This 
fluctuating birefringence causes the direction of polarization to fluctuate in phase with the 
ellipticity. The fluctuations of the carrier density have not been given here, since they are 
difficult to observe experimentally. They would be strongly correlated with the oscillating 
field terms, being 90 degrees out of phase with respect to the fluctuations of the field. 

Because of this correlation effect, fluctuations in the direction of polarization P 2 are 
always much stronger than fluctuations in the ellipticity P3 if a is greater than 1. Indeed, 
the values of 2 to 6 given for a in the literature suggest a difference of almost an order of 
magnitude. 

Another typical feature of the polarization fluctuations is that they approach infinity 
close to threshold. This is an indicator that the light emitted very close to threshold is still 
lamp like. A non vanishing amount of laser emission is necessary to overcome the noise 
effects and to stabilize both the intensity and the polarization. 

Although three unknown parameters, r, p and 9, enter into the model, the polarization 
noise terms calculated here are defined by only two parameters, namely p + 9 and r + p — 9. 
To fully separate the effects of spin relaxation and of anisotropy, one additional parameter 
is needed. This additional parameter may be found by taking a closer look at the nearly 
degenerate oscillation frequencies in the intensity and in the polarization. 



D. Perturbation theory for the frequency of relaxation oscillations 

The relaxation oscillations appear equally in the total intensity and in the polarization 
fluctuations. This is a direct result of the two density model: ellipticity fluctuations and 
the associated fluctuations in the polarization direction are a result of uncorrelated intensity 
fluctuations in the two circular polarizations. 
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However, anisotropies couple the two subsystems and induce slight changes in the fre- 
quency. By calculating the difference between the relaxation oscillations in intensity and in 
polarization, further information on the anisotropies can be obtained. 

The perturbative correction to the Eigenvalues may be obtained by calculating the matrix 
elements M i: j between the approximate Eigenvectors and hj. The correction to the 
Eigenvalue \ caused by a weak coupling to the Eigenstates of Xj is given by MijMj^/ '(Aj — 
Xj). Calculated to second order in 7/V, the Eigenvalue corrections are 



5A 2 



M li2 M 2il 
Ai — A2 
.7 2 x 2 



8u 
-5X 1 



(24) 

(25) 
(26) 



for the total photon number and 



5X 3 = 

.1/ ,..-,.!/;,..; M 4 , 3 M 3i 4 

OM — — ; r — — 

A4 — A3 



A4 — A5 

7 2 

—i — (x + r — p 
8u 



,2 , .7W + I) 



+ 1- 



2v 



5X?, = —5Xa 



(27) 
(28) 

(29) 
(30) 



for the polarization variables. The difference in frequency between the intensity oscillations 
v n and the polarization oscillations u P is 



Up = 



11 
8u 



(x + r — p+oy 



x 



49' 



,a 2 + 1 



a" 



(31) 



Note that only a frequency anisotropy will make up larger than u n . Therefore, if the po- 
larization noise oscillates faster than the intensity noise, this is an indicator of a strong 
frequency anisotropy. 



V. DISCUSSION OF THE RESULTS AND EXPERIMENTAL POSSIBILITIES 

A. Determination of timescales and anisotropies from the fluctuations of the laser 

light 

The equations given above show the wealth of information that can be obtained from 
measurements of the fluctuations in the laser light emitted from a VCSEL during cw oper- 
ation at a stable linear polarization. 

Intensity noise is given by the relaxation rate 'jx/2, the oscillation frequency u and the 
relative magnitude of Ajxix — 1). Experimental determination of these three quantities is 
equivalent to a measurement of the three timescales k(1 + 1), 7 and w(l+g). 

The ellipticity noise (P 3 (t)Ps(t + r)) differs from intensity noise in the relaxation rate, 
which is 7(2; + r + p — 9)/2 instead of jx/2, and in the relaxation frequency, although the 
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latter is only slightly different. These timescales also show up in the correlation of ellipticity 
and polarization direction fluctuations, (P 3 (t)P 2 (t + r)), which is just a times the ellipticity 
fluctuations. An observation of this correlation can provide not only strong evidence in 
support of the split density model, but is also a direct measurement of a. Finally, the 
fluctuations of polarization direction, (P 2 (t)P 2 (£ + t)), include not only noise correlated to 
ellipticity fluctuations, but also additional noise with a relaxation rate of 7(p + 9). From 
the two relaxation rates and the difference in oscillation frequency, the parameters r, p and 
9 can be calculated, thereby allowing an experimental determination of the spin relaxation 
rate, of the gain-loss aniosotropy and of the frequency anisotropy. 



B. A note on polarization stability 

The polarization fluctuations predicted by the split density model are largely fluctuations 
of the direction of linear polarization. Fluctuations of the ellipticity are smaller by at least 
a factor of a 2 . For very small anisotropies, the main contribution to the fluctuations will 
be from the weak relaxation rate of 7(p + 9), implying an even greater discrepancy between 
fluctuations in ellipticity and in polarization direction. 

To estimate the relative anisotropies necessary to stabilize polarization, it is therefore 
most appropriate to require the condition p + 9 ^> A(l + a 2 )/(x — 1) to be fulfilled. For 
reasonable estimates, p and 9 should then be at least of order 1. This would necessitate 
a relative gain-loss anisotropy g — I of at least 1CT 3 or a frequency anisotropy of at least 
1 GHz. 

However, these values depend critically on the timescales of the laser process, all of which 
enter into the magnitude of the fluctuations given by A. 



C. Experimental possibilities 

To determine the timescales and the anisotropies of quantum well VCSELs as described 
above, time resolved measurements of the polarization fluctuations at a resolution of at 
least picoseconds are required. Further, polarization filters for both linear and circular 
polarization are needed. 

Ideally, P 2 and P3 could be measured directly by separating the laser beam using a bire- 
fringent material and measuring the intensity difference between the two parts. However, it 
is also possible to measure the fluctuations by inserting a filter and measuring the fluctua- 
tions in n s (l + P 2 / 3 )/2. Since the fluctuations in n are not correlated with the fluctuations 
in polarization, the resulting relative noise is just the avarage of total intensity noise and 
polarization noise. E.g. if one measures the fluctuations in the intensity of right circular 
polarized light one obtains 

</+(l)/+(l + T» = U WW + T)) +m)pt(t + T)) \ 
1+ ^ \ U s ) 

A measurement of the correlation of ellipticity and polarization direction fluctuations is 
very difficult, since it requires a splitting of the laser beam without causing undesireable 
polarizations. Probably the best approach would be to reflect only a small fraction of the 
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beam at a right angle, using this light for the measurement of linear polarization direction 
and using the nearly unaltered beam for the measurement of ellipticity fluctuations. Another 
possibility is to keep track of all the polarization properties of the optical devices used, 
sorting out the contributions to the measured noise afterwards. The Stokes parameters 
are quite convenient for this task, since all polarization effects can be described by matrix 
multiplication. 

In this manner, all timescales and anisotropies can be determined at a fixed injection 
current. If such an experiment works, it is also possible to measure changes in these prop- 
erties as injection current is increased. A comparison with the temperature dependence at 
constant injection current is also possible. This should reveal a lot of the material properties 
of the semiconductor and help to clear up some of the open questions regarding polarization 
switching. 



VI. CONCLUSIONS AND OUTLOOK 

The analysis presented here shows how all the parameters of the two density model 
for a VCSEL with stable linear polarization can be obtained from measurements of the 
polarization fluctuations. Time resolved measurements of the polarization fluctuations in 
VCSEL light can therefore unambiguously resolve open questions, such as whether temper- 
ature effects cause polarization switching or whether birefringence or gain-loss anisotropies 
are responsible for polarization stability. 

This is extremely important, because experimental results on polarization and intensity 
as a function of the injection current can only be interpreted when the correct mechanism of 
polarization stability has been identified. For example, the authors of assume that only a 
frequency anisotropy contributes to polarization stability without considering the possibility 
that a gain-loss anisotropy might have an effect as well. Unverified assumptions are also the 
subject of the criticism voiced in || about the explanations given in ||. Further experiments 
such as the ones proposed here are absolutely necessary to avoid misleading interpretations. 

The polarization fluctuations show features typical for the split density model, which 
can be used as a test criterium for whether the model is valid in a given device or not. 
Since the split density represents the effect of quantum well confinement on the polarization 
properties of VCSELs, this effectively tests the quantum well structure in the active region. 
The quantity which depends strongly on the size of the quantum wells is the spin relaxation 
rate. Experimental results on this rate are also of interest in connection with calculations 



II4II6I and luminesence experiments [|13j carried out to investigate spin flip scattering in 
quantum wells. 

The rate equations presented here are formulated in a very general way and may also be 
applied to cases with more excotic anisotropies and timescales, such as discussed in || and 
||, which show switching and/or include a magnetic field. In all these cases, the investigation 
of noise adds additional predictions for experiment to the results and thereby increases our 
understanding of the physics involved in the polarization properties of VCSELs. 
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FIGURES 



FIG. 1. The arrows on the Poincare spheres illustrate the dynamical effect of anisotropics on 
the normalized Stokes vector P. The sphere on the left shows how the frequency anisotropy causes 
P to rotate around the axis defined by the anisotropy. The sphere on the right shows how the 
gain-loss anisotropy pulls P towards one of the poles defined by the vector of anisotropy. 

FIG. 2. Contour plot of the Gaussian Distribution corresponding to the polarization fluctua- 
tions atr = 0forx = 2,a = 2,r = 2,p = 2 and = 2. This choice of parameters clearly shows 
the correlation between polarization direction and ellipticity. For A = 0.01, the fluctuations of P2 
correspond to deviation of approximately ±5 degrees in the direction of polarization. 

FIG. 3. Time dependence of the fluctuations (Pi(t)Pj(t + r)) for x = 2, a = 2,r = 2, p = 2 
and 6 = 2. As explained in the text, the fluctuation of polarization direction = (2,2) is 

largest, while the correlation of polarization direction and ellipticity = (2,3) is exactly equal 
in magnitude to two times the fluctuations in ellipticity = (3,3). Time is given in units of 
7 -1 , which is typically 100 ps to 1 ns. The variable A is typically around 1/100. 
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